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Grain boundaries are topological defects that often have a disordered character. Disorder implies 
that understanding general trends is more important than accurate investigations of individual grain 
boundaries. Here we present trends in the grain boundaries of graphene. We use density-functional 
tight-binding method to calculate trends in energy, atomic structure (polygon composition), chemical 
reactivity (dangling bond density), corrugation heights (inflection angles), and dynamical properties 
(vibrations), as a function of lattice orientation mismatch. The observed trends and their mutual 
interrelations are plausibly explained by structure, and supported by past experiments. 

PACS numbers: 61.72.Mm,61.48.Gh,73.22.Pr,71.15.Nc 



I. INTRODUCTION 

Real graphene has always defects, edges, point de- 
fects, chemical impurities — and grain boundaries. Grain 
boundaries (GBs) are topological defects, trails of dis- 
order that separate two pieces of pristine hexagonal 
graphene sheets. Th^ reside practically in any graphitic 
material, graphitj^^, soo lP^^^, single- and multi-layer 
graphend^, fuUereneJ^EHl^ carbon nanotube d^^ * ^^!, with 
varying abundance. In graphene fabrication, for instance, 
chemically synthesized samples tend to have more GBs 
than mechanically cleaved samples^^. For electron mobil- 
ity high purity is important and GBs are best avoided, 
while other properties may tolerate defects better. We 
refer here to strictly two-dimensional GBs (albeit not 
necessarily planar) which should not be confused with 
genuinely three-dimensional GBs in graphite. 

But not are defects, impurities and GBs always some- 
thing bad; they are interesting also in their own right — 
even usefuP^EH. After all, being extended, GBs mod- 
ify graphene more than point defects. They affect 
graphene's m agnet ic, electronic, structural, and mechan- 
ical propertie d l ^^l 

Considering the relevance and prevalence of graphene 
GBs, they ought to deserve more attention. Most re- 
lated studies are on GBs in graphite and moire pat- 
terns thereirP^^it^, or on intramolecular junctions in 
carbon nanotubedi^ttSli qbs in si ngle- or few-layer 
graphene, whether experimentaP^'^or theoreticaP^, are 
still scarce. Often graphitic grai n bou ndaries are ide- 
alized by pentagon-heptagon pair d^^*^^ , where the pair 
distances determine lattice mismatch. While this ideal 
model works in certain occasion d^^ * ^^ !, it does not work 
for GBs that are rough, corrugated, have ridgefl^, and — 
most importantly — are decorated by dangling bonds or 
reactive sites alik^. 

In this paper we investigate graphene GBs beyond sim- 
ple models. We explore computationally the trends in 
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structural properties, chemical reactivity and vibrational 
properties for an ensemble of GBs that are not ideal, 
but have certain roughness. What this roughness really 
means will be clarified in Sec. III. 



II. COMPUTATIONAL METHOD 

The electronic structure method we use to simulate 
the GBs is density-functional tight-binding (DFTB jSlMl^ 
and the hotbit codd^. DFTB models well the covalent 
bonding in carborpSH2n!^ ^nd suits fine for our simula- 
tions that concentrates on trends. We leave the details 
of the method to Ref. US] and only mention — this will 
be used later — that the total energy in DFTB, includ- 
ing the band-structure energy, the Coulomb energy, and 
the short-range repulsion, can be expressed as a sum of 
atomic binding contributions, 

N 

£^DFTB = ^ Ci- (1) 
i=l 

Here N is the atom count and atom i's contribution to 
cohesion energy; in pristine graphene = egj. — —9.6 eV 
for all atoms (DFTB has some overbinding) . The quan- 
tity Ci allows us to calculate energies in a local fashion: 
for a group of carbon atoms in a given zone S (set of 
atoms), the sum 

Ai? = ^(Q-egr) (2) 
ies 

measures how much the energy of zone S is larger than 
the energy of equal-sized piece of graphene. In addition, 
since egr = —9.6 eV comes from 3 equivalent bonds, the 
value ei — 2/3 • egj. « —6.4 eV infers a dangling bond 
for atom i. Hence, albeit less accurate than density- 
functional theory, DFTB enables straightforward anal- 
ysis of local quantities, in addition to fast calculations 
and extensive sample statistics. 
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FIG. 1: (color online) Grain boundary sample, one out of 48, 
made by merging two graphene edges with the same chirality, 
but different orientation and random translation along the 
boundary; our samples are hence characterized by chirality. 
In this figure chiral indices are (11, 6), chiral angle 9 — 20.4°. 
Dashed vertical lines stand for periodicity across the horizon- 
tal direction. 



III. SAMPLE CONSTRUCTION 

To investigate trends, we constructed an ensemble of 
48 graphene GBs using a procedure we now describe, (i) 
Cut two ribbons out of perfect graphene in a given chiral- 
ity direction, with random offset. Chirality index (n, m) 
means we cut graphene along the vector C — nai +ma2, 
where an are the primitive lattice vectors. The cut direc- 
tion can be expressed, equivalently, vy the chiral angle 
9 = tan^^ [-\/3to/("i + 2n)]. (ii) Passivate ribbons' other 
edges with hydrogen, leaving carbon atoms on the other 
edges bare, (iii) Thermalize the ribbons, still separated, 
to 1500 K using Langevin molecular dynamics (MD). (iv) 
Perform reflection operation on the other ribbon with re- 
spect to a plane along ribbon axis and perpendicular to 
ribbon plane, (v) Perform random translations along the 
ribbon direction, (vi) Merge the two ribbons' bare car- 
bon edges with MD at 1500 K. (vii) Cool the system 
to room temperature within ^ 1 ps. (vii) Optimize the 
geometry using the FIRE methocPI^, and finally obtain 
structures like the one in Fig[T] The total process takes 
about 3 ps. 

Our ensemble of 48 GB samples are chosen to get rep- 
resentatives for all chiral angles from = 0° to 6* = 30°. 
Our limit for the width in cutting the ribbons is 10 A. 
Ribbons have to be wide enough (> 9 A) for faithful de- 
scription of a semi-infinite graphene, and, by limiting the 
total number of atoms below 170, also the periodicity of 
GB becomes limited below 40 A. (Edges with 6* ^ 0° and 
9 ^ 30°, that would have short periodicity, are scarce.) 
Our procedure to make GBs is neither the only nor the 
best one, but is suitable for trend-hunting. 

Sure enough, it's a downside that the procedure ex- 
cludes GBs from ribbons with different edge chirality. 



such as merged zigzag and armchair edges. Different 
chirality, unfortunately, would mean different periodic- 
ity and practical problems. To have a single number, 9, 
to identify GBs is vital for finding the trends we con- 
centrate on. For GBs identified by two numbers, 9i and 
92, trend- hunting would be harder; those GBs can be 
investigated afterwards, using the insights we learn here. 
Further, while the temperature 1500 K, motivated partly 
by experiments^, allows rearrangements in merging, the 
heat of fusion renders initial temperatures irrelevant. We 
use canonical MD in the merging to allow GB formation 
process to he steered by intrinsic driving forces and ener- 
getics. 

The time scale of the process is limited, as usual, in 
part by computational constraints. However, prolong- 
ing the process would cause no fundamental structural 
changes as the merging process itself is nearly instan- 
taneous; there is no diffusion after merging and satura- 
tion of the dangling bonds, there is only annealing of the 
worst local defects and stress-release that causes buckling 
(Sec VI). The usual time scales in graphene growth are 



hence not relevant to our GB formation process. Only in- 
trinsic roughness remains in the GB samples — roughness 
related to chiral angle and randomness in cutting offsets 
and translations. 

Usually all atoms in ribbons were also bound to GBs, 
but in few structures (near armchair chirality) the merg- 
ing process squeezed dimers out, and either inflected 
them out of plane or detached them completely. Since 
ribbons found these own ways to optimize geometries — 
ways not anticipated beforehand — it suggests that the 
design of our construction process is not dominant. 



IV. PRELUDE: ENERGIES FOR FREE 
GRAPHENE EDGES 

Before going into GBs themselves, let us look at the 
free bare edges of graphene. Fig[2]shows the edge energies 
as a function of the chiral angle, 9 = Q° meaning zigzag 
and 9 = 30° armchair edges. To ignore the effect of 
hydrogen passivated edge, carbons bound to hydrogen 
are neglected in edge energy calculation, 



^edgc 



L ^ 

iGcdgc 



(e, -egr). 



(3) 



The edge energy between zigzag and armchair varies 
linearly; fluctuations in energy are due to random offset 
in the cut, occasionally producing pentagons. The edge 
energy in zigzag is high due to strong and unhappy dan- 
gling bonds; in armchair the dangling bonds are weak- 
ened bythe formation of triple bonds in the armrest 
partPEl. 



It was recently predicted theoreticalljPZl and later con- 
firmed experimentallji^ that the zigzag edge is actually 
metastable, and prefers reconstruction into pentagons 
and heptagons at the edge, forming a so-called reczag 
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FIG. 2: (color online) Energies of the free edges of graphene as 
a function of chiral angle; zigzag edge has = 0° and armchair 
6 = 30° . Free zigzag edge segments are metastable and prefer 
local reconstruction (two adjacent hexagons reconstruct into 
a pentagon and a heptagon); these so-called reczag edges are 
shown with gray connected symbols. The reczag segments 
are, however, important only for free edges and irrelevant for 
grain boundaries studied here. 

edge. As it turned out, reczag is energetically even bet- 
ter than armchair. Hence, if we reconstruct the zigzag 
segments in edges with small 0, we usually lower the edge 
energy, as seen in Fig [2] where reconstructions are added 
by hand. Reczag edges are, however, irrelevant for our 
GBs where edges are not free, and are ignored because we 
want the dangling bonds to spontaneously find contact 
from the other merging edge. The edge energies were in- 
vestigated and presented here for comparison with DFT 
calculations^^. The accuracy in edge energies is better 
than 10 %, and we expect same accuracy in GB energet- 
ics. 



V. TRENDS IN ENERGY AND STRUCTURE 

The GB energy per unit length is 

'^GB = XI ~ ^sr) : (4) 

ieGB 

measuring how much GB costs energy relative to the 
same number of carbon atoms in graphene. The energies 
for the whole ensemble of GBs are shown in Fig[3j as a 
function of the chiral angle. The value Egb — 0.0 eV/A 
appears only with 9 = 0° and 9 = 30° — meaning pristine 
graphene. 

Most GB energies are less than or equal to the edge 
energies of the free graphene edges, albeit with variation. 
This means that GBs regain, on average, other ribbon's 
edge energy during the merging, and the energy of fu- 
sion ranges from 1 ... 2 eV/A; on average half of the free 



FIG. 3: (color online) Trends in energy density for graphene 
grain boundaries as a function of the chiral angle. The zero- 
energy means pristine graphene, and is achieved only for 
zigzag (6 = 0°) and armchair (9 = 30°) chiral angles. 

edges' dangling bonds get passivated. The picture is not 
this simple, however, as MD simulation creates different 
polygons that cause strain. The randomness of the poly- 
gon formation is manifested by the energy variations in 
Fig]3] Compared to ideal GBs with pentagons and hep- 
tagons only, our ensemble reveals the full complexity that 
rough GBs have. 

GB energies are the highest around 6* ^ 15°, which can 
be understood from structural analysis, shown in Fig[4^ 
as polygon distribution within the GB zone. What we 
count into GB zone are all the polygons that were not 
part of the ribbons prior to merging. Polygons larger 
than nonagons, which appear more like vacancies instead 
of polygons, are omitted here. It is around 6 ~ 15° where 
the abundance of hexagons is at minimum and GBs are 
invaded by other polygons. The abundance of pentagons 
and heptagons is as high as the abundance of hexagons, 
but also squares, octagons and nonagons are found. Close 
to 6* ~ 0° and 9 ~ 30° hexagons prevail. From this we 
may conclude that around 9 ~ 15° the edge geometries 
have the largest mismatch, resulting in various polygons, 
consequent strains, and high energy. 

To understand how much different polygons cost en- 
ergy, we used the quantity 

En = ^ - Cgr), (5) 

measuring how much atoms, on average, cost more in n- 
gon relative to atoms in pristine graphene; this quantity 
is for illustration only and considers polygons as sepa- 
rate items — the sum of e„'s for given GB is not the to- 
tal energy. Fig|4|3 shows e^'s for the polygons within 
GB zones as a function of the area of the polygon (de- 
termined by triangulation) . For simple geometrical rea- 
sons for small polygons the area distribution is narrow; 
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FIG. 4: (color online) a) Relative number of polygons in 
grain boundaries as a function of chiral angle. The largest 
lattice mismatch, near chiral angle ~ 15°, is manifested by 
the appearance of mainly pentagons and heptagons, but also 
squares, octagons and nonagons. In b) each dot represents one 
polygon in all the grain boundary polygons in all of our 48 
samples. Abscissa shows the area of that polygon (determined 
by triangulation) and ordinate shows the polygon energy cost 
per atom, relative to an atom in pristine graphene. 



VI. TRENDS IN BUCKLING (INFLECTION 
ANGLES) 

In the generation process GBs are free to deform. The 
optimized GBs will typically end up having inflection an- 
gles, as the inset in FigjS] illustrates. The data points in 
FigjS] show the inflection angles for the GB samples as 
a function of chirality. Flat GBs occur with 9 — 0° and 
6 — 30°, that is, with pristine graphene alone. 

The notable trend in this scattered plot is the system- 
atically small inflection angles, meaning flat GBs, around 
6 ~ 30°; near 9^0° inflection angles are more scat- 
tered meaning geometries that vary from flat to sharply 
kinked GBs. This trend can be understood by edge pro- 
files: the stronger dangling bonds at zigzag edge, when 
brought into contact with another edge, can induce larger 
distortions than inert armchair edge (dangling bonds at 
armchair are partly quenched by triple-bond formation) . 
The steps in edge profiles near ~ 0°, moreover, are 
2.1 A, whereas the steps in edge profiles near 9 ^ 30° are 
only 1.2 A. This means that, to make bonds, the edge 
atoms near 9^0° need more pulling than edge atoms 
near 9 ~ 30°, causing the buckling. 

It is clear that the numbers in FigjS] have no direct 
relevance as such, because our geometrically optimized 
GB samples are in vacuum, and measure only ^ 20 A 
across the GB. We argue, however, that the inflection 
angle measures how much given GB would buckle in an 
experiment — large inflection angle meaning tendency to 
stick out (sticking out would be bound by geometric con- 
straints and by surface adhesion). For example, a GB 
with 9 ^ 30° on a support will always remain smooth, 
whereas a GB with ~ 0° will either remain smooth or 
buckle, to be seen as a ridge or a mountain range in scan- 
ning probe experiments.lSl Indeed, this buckling tendency 
was seen in an experiment by Gervenka at al. , measuring 
corrugation heights up to 3 A within GB region^. Small 
corrugation heights and inflection angles are natural in 
samples on flat substrates, but also large inflection angles 
are realistic, for example in soot particles.!^ 



VII. TRENDS IN CHEMICAL PROPERTIES 



large polygons have more freedom to change their shape. 
For small polygons the distributions in e„, on the con- 
trary, are wider; this is partly due to smaller n in Eq.([5]). 
Clearly, the cost e„ of any polygon depends on its envi- 
ronment, just as it also depends for hexagons, for which 
eg = ... 1 eV. But note that Figjijs already contains 
the effect of the polygon environment and all potential 
cross-correlations (such as pentagons often neighboring 
heptagons). It would be interesting to investigate poly- 
gon statistics also from transmission electron microscopy, 
now that aberration-corrected measurements can achieve 
atom accuracy^. 



Now we pose the general question: how reactive are 
GBs? We approach this question by examining hydro- 
gen adsorption energies and dangling bonds (DB). To 
this end, we first develop a connection between hydrogen 
adsorption energy and the electronic structure given by 
DFTB. 

Figj6|D shows the hydrogen adsorption energies for car- 
bon atoms in 6 representative GB samples, from total 
350 hydrogen adsorption calculations. It appears that 
adsorption to given carbon atom i is strong if atom's 
cohesion — decreases. Carbon atoms part of regular 
hexagons have adsorption around —2 eV (® in Fig[6^) 
but adsorption increases if atom is surrounded by other 
polygons and bonding angles deviate from 120° (@, (D, 
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FIG. 5: (color online) Inflection angles after optimizations for 
all grain boundary samples; 0° means planar structure (oc- 
curring only with pristine graphene, 9 — 0° and 6 — 30°). 
The error bars remind of technical ambiguities in angle de- 
termination. The main trend here is that zigzag edges cause 
both large and small inflection, whereas armchair edges cause 
only small inflections. Inset: inflection angle illustrated. 
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and ® in Fig|6p). The common denominator in these 
examples is that the change in hydrogen adsorption en- 
ergy is caused by the strain in bond angles and in bond 
lengths. 

The strongest adsorptions around —5 ... — 6.5 eV, in 
turn, are caused by dangling bonds that have ei > — 7 eV 
(® in Figj6]); such strong adsorption never occurs with 
three-coordinated atoms. (We gave the argument about 
DB energetics already in Sec|llj) We note that DFTB hy- 
drogen adsorption energy to zigzag (5.8 eV), for example, 
agrees reasonably with DFT energy (5.4 eV)^^. There- 
fore, we can characterize the reactivity directly by the 
DFTB electronic structure using the quantities e^, with- 
out any adsorption calculations. Surely, dangling bonds 
can be identified from geometry by defining criteria for 
coordination numbers, but this approach is prone to er- 
rors, especially in disordered regions that GBs are. 

Using ei > —7 eV as a criterion for a dangling bond, we 
then analyzed the reactivity for the whole GB ensemble. 
Fig|6]: shows the average number of dangling bonds per 
unit length in a GB with a given chiral angle. The highest 
density, one DB per ~ 14 A, occurs with 9 ^ 20°. The 
lowest density, one DB per ~ 30 A, occurs with 6 ^ 15°; 
pristine graphene with no DBs becomes probable with 
6* - 0° and e - 30°. 

By analyzing the histogram in FigjOj: another way, 
30 % of the samples have no DBs, 54 % of the samples 
have one DB every 7.5 ... 20 A, and 16 % of the sam- 
ples have one DB every 20 ... 50 A. These numbers agree 
with a recent scanning tunneling microscopy (STM) ex- 
periment. Namely, dangling bonds are highly localized 
states just below the Fermi-level, and hence seen as a 
bump in constant-current STM with low bias. The STM 
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FIG. 6: (color online) Characterizing reactivity in grain 
boundaries, a) Grain boundary sample with (9, 5) chirality, 
including various polygons, b) Characterizing chemical reac- 
tivity from the DFTB electronic structure directly: we cal- 
culated hydrogen adsorption energy for selected 350 atoms as 
adsorption sites and plot for each atom adsorption energy and 
€i. Given the correlation between the adsorption energy and 
Ei, we access reactivity directly from the electronic structure, 
without additional calculations. The adsorption site numbers 
refer to panel a); atom CD, having three ~ 120° angles has a 
small adsorption energy, atom ® has a dangling bond, and 
other atoms have strained environment, c) Using the correla- 
tion established in panels a) and b), we calculate the averaged 
density of dangling bonds per unit length along grain bound- 
ary for all samples, as a function of 9. 
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FIG. 7: (color online) Projected vibrational density of states 
(PVDOS) as a function of the wavenumber and the chiral an- 
gle. PVDOS for each sample (given 8) was calculated by solv- 
ing vibrational spectrum from the dynamical matrix, and pro- 
jecting the vibration eigenmodes on atoms within the grain 
boundary; the contour plot is then gathered and smoothed 
from all the 48 samples. Intensity increases from blue to red. 



some 100 cm~^ in energy around 9 ^ 15°. This can be 
explained by structure. As discussed in Sec|Vj around 
9 ^ 15° hexagons are at relative minimum, implying a 
non-uniform and less rigid structure, and causing floppi- 
ness in high-energy modes. The high-energy modes have 
bond stretching between nearest-neighbors and are hence 
sensitive to local structural changes. Low-energy modes, 
again, are more collective and hence insensitive to local 
changes in structure, as long as they remain somehow 
"graphitic" . If FigjT] and Fig|4] are compared carefully, 
one finds that the local intensity maxima of the G-band 
occur precisely for 9 with local maxima in the abundance 
of hexagons. The intensity increase around ~ 750 cm~^ 
and 6* ~ 15° is mainly due to renormalization of PVDOS. 

The observations above are qualitatively supported by 
earlier experiments. Using Raman spectroscopy, a G- 
band shift from 1600 cm^^ to 1510 cm^^ was reported by 
Ferrari and Robertson in Rcf. 35; the shift was identified 
to be a consequence of structural change from nanocrys- 
talline graphite to amorphous phase. While still being 
mainly planar (apart from inflection), GB zone around 
9 ~ 15°, such as Figj6^ with 9 = 20.6°, indeed appears 
amorphous. 



images of Cervenka et al. in Ref. H] show periodic ap- 
pearance of sharp peaks, with 5 ... 20 A periodicity for 
53 % of the samples and with > 20 A periodicity for 
47 % of the samples. Even if these periodicities should 
be caused by adsorbed impurities and not from bonds 
that dangle, it is still likely to be result from reactive 
sites within GBs — and answers the original question we 
posed in this section. The agreement with experiment 
gives confidence in the objectivity of the GB construc- 
tion process. Understanding the trends in defects and 
reactivity will hopefully help in the design of functional- 
ized graphene compounds^. 

VIII. TRENDS IN VIBRATIONAL 
PROPERTIES 

The structure of a GB, given its constituent polygons, 
affects directly on its vibrational spectrum, and gives ex- 
perimentally complementary information. 

FigjT] shows the projected vibrational density of states 
(PVDOS) for the GBs as a function of chiral angle and 
wave number. PVDOS was calculated by solving vibra- 
tions for the whole GB, by projecting the eigenmodes to 
atoms within GB area, and by renormalizing the modes — 
for all the 48 GB samples. Hence FigjTjshows a lot of data 
in a compact form. 

Spectra show three main bands, two at low ener- 
gies ~ 750 cm~^ and ~ 500 cm^^, and one — the so- 
called G-band in Raman spectroscopy — at high energy 
~ 1600 cm~^. The band at ~ 500 cm~^ is steady across 
all 0's and can not be used for structural identification. 
The G-band, in turn, loses intensity and comes down 



IX. CONCLUDING REMARKS 

Some defects, like singular Stone-Wales defect or 
reczag edge of graphene, have their own name because 
they are well defined. Grain boundaries, on the contrary, 
are like snowflakes — there is no flake like another, but it's 
enough to know that they are roughly hexagonal, flat, 
small and cold. For the same reason it's enough to know 
how grain boundaries usually look and feel. Knowing 
trends is valuable. 

We investigated GBs from different viewpoints, discov- 
ering trends with complementary information, measur- 
able also experimentally. We investigated (i) geometry 
(polygons and inflection angles) measurable with trans- 
mission electron microscope or atomic probe microscope; 
(ii) energy, that is manifested in geometry and thereby 
measurable; (iii) reactivity, measurable with adsorption 
experiments; (iv) vibrational properties, measurable with 
Raman spectroscopy. We are confident that the trends 
are genuine, part because of agreements with earlier ex- 
periments, part because the trends all make intuitive 
sense: energy, inflection angles, reactivity and vibration 
trends make sense given the structure, and the structure 
trends make sense given the graphene edge mismatch. 
We hope the trends help to explain experiments — and 
also help simulations and experiments to design graphene 
structures for given functions. 
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